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Abstract 

We consider the adiabatic evolution of Kramers degenerate pairs of spin states in a half-integer 
spin molecular magnet as the molecule is slowly rotated. To reveal the full details the of the 
quantum evolution, wc use Majorana's paramctrisation of a general state in the 2j + 1 dimensional 
Hilbert space in terms of 2j Majorana points. We show that the intricate motion of the Majorana 
points may be described by a classical hamiltonian which is of the same form, but of quite differ- 
ent origin, as that which appears in the spin-coherent-state path integral. As an illustration we 
consider molecular magnets of the j = 9/2 Mn4 family and compute the frequency with which the 
magnetization varies. This frequency is generally different from the frequency of the rotation. 

PACS numbers: 75.50.Xx, 03.65.Vf, 03.65.Aa 
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I. INTRODUCTION 



The controlled and measurable manipulation of quantum spin states is the basis of NMR 
with its abundance of practical applications. Another system in which spins may be ma- 
nipulated and measured is provided by molecular magnets [1]. These materials consist of 
large organic molecules containing clusters of transition elements such as Mn or Fe, whose 
unpaired electron spins lock together to form a single large (j ~ 10) spin. The best known 
examples are Mnl2, Fe8, and Mn4, for which j — 10, 10 and 9/2 respectively. (The nomen- 
clature focusses on the magnetic atoms. Mnl2, for example, is shorthand for the polyacetate 
Mni20i2(CH3COO)i6-(H20)4.) The organic ligands surrounding the magnetic cluster serve 
to reduce the interaction between clusters so that each molecular spin behaves almost in- 
dependently — yet, when the molecules have a common orientation and all experience the 
same external forces, the total magnetic moment is the sum of many identical molecular 
moments and so the spin direction may be observed. The energy of the spin J depends on 
its orientation with respect to the host molecule, and one may think of the tip of the spin 
arrow as moving on a sphere through a potential landscape possessing hills, valleys and sad- 
dlepoints. This potential landscape can be manipulated by applying an external magnetic 
field. In Fe8, for example, the rate of quantum spin tunneling between eigenstates located 
in separate potential wells has been measured via the Landau-Zener crossing of the levels as 
the external field is varied [2]. 

An alternative route to manipulating the spin is simply to rotate the system. In most 
cases the magnetization will merely co-rotate with the molecule, but when the spin ground- 
state is degenerate we have more interesting possibilities. If the rotation is slow compared 
to the separation of the degenerate ground state from the first excited state, the evolution 
will be described by non-abelian Berry transport [3, 4]. Now a doubly-degenerate ground 
state is guaranteed by Kramers' theorem whenever the total spin is a half integer, and the 
Hamiltonian H{J) governing the large spin's interaction with the host molecule is invariant 
under time reversal [5]. The aim of this paper is to explore how these Kramers-degerate 
pairs evolve as a general anisotropic molecule is rotated. 

When dealing with the 2j + 1 dimensional Hilbert space of a large spin, the expectation 
(iplJlip) reveals only a small part of the quantum information encoded in the spin state {ip). 
A useful tool [6] for visualizing the complete quantum state is provided by the location of 
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the zeros on the unit sphere of the spin-coherent-state wavefunction ip{z). These zeros are 
the antipodes of "Majorana points" that can be thought of the directions of a set of 2j 
spin-l/2's that compose the spin-j. 

In section II we will review Majorana's parametrisation [7] of a general spin-j state 
and explain its connection with spin-coherent state wavef unctions. In section III we will 
derive expressions for the non-abelian Berry connection one-form in terms of the Majorana 
parametrisation, and show how they simplify when we restrict to Kramers degenerate pairs 
of states. In sections IV and V we introduce an alternative way of thinking about the 
adiabatic Berry transport for Kramers pairs, and how the evolution of the state is described 
by a classical hamiltonian dynamical system. In the last section VI we show how even simple 
anisotropics can give rise to intricate state evolutions as a molecule is rotated, and that this 
intricacy may have experimentally observable consequences. In particular we will see that a 
steady period T rotation of the molecule can give rise to motions of the spin direction that 
are again periodic, but with a period that typically differs from T. 

II. MAJORANA PARAMETRISATION AND SPIN COHERENT STATES 

It is a familiar fact that when a collection of 2j spin-1/2 particles are combined into 
a single quantum system the spin-j representation of SU(2) resides in the space of totally 
symmetric tensors. Elements of a symmetric tensor space Sym [V"^^] are sums of the form 

A''''-'^ei,Qei,Q---Qei^, (1) 

where the are the basis vectors for V, and the symmetric tensor product "0" is defined 
to commutative and distributive: 

a0b = b0a, a (Ab + /xc) = A (a b) + (a c). (2) 

In [7] Majorana showed that an arbitrary spin-j state ^ e Sym[\/^y2''] can be decomposed 

as 

* = ^,,0---0V'i,,, (3) 

where the 

t/ji = aiei+ /3i 62 (4) 
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are a set of 2j spin-1/2 spinors with ei = |t), 62 = This result is perhaps surprising. 
Most elements of a general tensor product space (whether symmetric or not) cannot be 
decomposed into a single product; they can only be expressed as a sums of products. Majo- 
rana's result was rediscovered by Schwinger [8] and by Rabi [9], and later reviewed by Bloch 
and Rabi [10]. Schwinger's paper is historically interesting as in his equation (9) he draws 
attention to a previously omitted term in the equations describing the evolution of a spin 
in a time-varying field. This term is now recognizable as the "Berry Phase." Majorana's 
original paper is similarly notable in that it contains an independent derivation of the non- 
adiabatic level-crossing probability that is traditionally attributed to Landau, Stuckleberg 
and Zener. 

Majorana's recipe for his decomposition is as follows: Given a spin-j state 

* = l^') = J]|j,m)(j,m|^') (5) 

m 

we construct a set of 2j + 1 coefficients 



cj.m = {j,m\^)^l m = -j,-j + l,...,j-l,j (6) 

and the degree 2j polynomial equation 

P^{z) = J2z'^--c^ = 0. (7) 

n=0 

If the 2j roots of this polynomial are denoted by Zi, then the Zi are related to the coefficients 
of the spin-1/2 factors by 

Zi = -Pi/ai. (8) 

We can extract the individual coefficients from their ratios by normalizing so that 
lo'iP + lAP = 1) but a choice of phase remains for each '0^. These phases can be swapped 
amongst the factors leaving only an overall phase to be matched with that of 

It may happen that one or more of the leading coefficients c„ are zero. The polynomial 
equation then has degree less than 2j and possesses fewer than 2j roots. In this case it 
is convenient to extend the complex plane to the Riemann sphere, and regard the missing 
roots as being at infinity. Their corresponding ctj are then zero. 



5 



To understand why the recipe works, we use the commuting property of the symmetric 
tensor factors to allow the shorthand 



ei • • • ei 62 • • • 62 = 6i+"'e^ (9) 

^ V ' ^ V ' 

j+m fcLctors j—m fcLctors 

To connect an expansion in terms of the e{'^"^e2~"^ with the more conventional expansion 

= J] |j,m)(j,m|*) (10) 

m 

we need to define an inner product on the symmetric tensor-product space. The one that is 
induced naturally from the spin- 1/2 product sets 

/oJ+"^o-'~"*loJ+"^'o-'~"^'\ — A (j + m) ! ( j — m) ! . . 

t=2 1*^1 / — "mm' ' 



Prom this we identify 



l^'' ^) = \l r^ ^Mr M^""^""- (12) 

^' [j + m.)\{j — my. 



This natural inner product arises from regarding the symmetric tensor products as being 
symmetrized conventional tensor products. For example, 

61 6i 62 = (6i 6i 62 6i 62 6i 62 6i 61) . (13) 

In the general case, the 3 is replaced by 

N. = M (14) 

^"^ (j+m)!(j-m)!' 

which is the number of distinct terms generated when we symmetrize 

6i, 6i, ■ ■ • © Bi^. = ^Vd) ® ^i.(2) ® ■ ■ ■ ® (15) 

by summing over the (2j)! elements of the permutation group &2j. The Njm distinct terms 
are mutually orthonormal with regard to the standard inner product on the unsymmetrized 
[^i/2]'^^'' , so computing the inner product from this conventional tensor- product side makes 

(6f "^6r"^|ef'"6r'") = {Njm)-%m = {N,m)-\ (16) 
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To confirm the identification, observe that J_ei = 62, J- 62 = 0, and so 



" [j + m)\[j — m)\ 



(j + m)!(j -m) 



■(j +m)ef 



(2j)! ^ (j+m-l)!(j-m + l)! ,. 



= A/0' + m)(i-m + l)|i,m-l) (17) 

as it should. 

Now let = ccjei + /3je2. Set Zi — —^i/ai, and expand out a decomposable state as 

* = • • • ^i,, = {\[ j (18) 
where a^-^ is the (j — m)-th elementary symmetric function of the —Zi. In other words 

n=0 i=l 

We therefore have 

ami^>=fn..V,_,,. /^'+'"^!;]j,-'")' . (20) 



or 

2i 



The common factor Hiii — (JjJ I^) docs not affect the location of the zeros, so we can 
replace the Oj-^ by Majorana's Cj_^ (which are (j, j|^)aj_^). Thus given the components 
{j,m\^) we can solve for the roots of P-ii,{z) — Yl^n=o^'^^~^^n — and construct the ctj, 
up to an overall phase as before. 

The Majorana polynomial P^^{z) for a given has a physical interpretation as the 
holomorphic spin-coherent-state wavcfunction representation of |^'). This wavefunction is 
obtained by first defining a family of spin coherent states 

\z) = exp{zJ+)\j, -j), {z\ = {j, -j\exp{zJ_) = \zy. (22) 

These states are not normalized, but have the advantage that the {z\ are holomorphic in 
the parameter z — i.e., they depend on z but not on z. The coherent-state wavefunction 



il){z) =^ {z\ip) corresponding to a state is then a holomorphic function. The label z can 
be identified as the complex sterographic co-ordinate of the spin direction on the two-sphere, 
with z = corresponding to the south pole (spin down) and z = oo to the north pole (spin 
up). (Strictly speaking |oo) is undefined, and one needs two families of coherent states with a 
transition function to capture every spin direction — see for example [11]. This technicality 
is not relevant for what follows however.) 

The inner product of two states and \x) may be evaluated in terms of their wavefunc- 
tions as 

,,, , 2j + l f dzAdz , , , , , 

Normalizable wavefunctions are therefore polynomials in z of degree less than or equal to 
2j. In particular 

{^\J, m) = J ^.^f/' z^^-\ (24) 

V -m)\{j +m)\ 

and the wavefunction of coherent state |C) is the polynomial 

{z\C) = {I + zCr^ . (25) 

When acting on the holomorphic wavefunctions then su(2) generators J3 and J± — Ji±i J2 
become 

J+ -z^— + 2jz 

J ^ — ^ 
dz^ 

>^3^^^-J- (26) 

It is straightforward to verify that with respect to the inner product (23) we have = J3 
and J± — J^:- 

Now from (6) and (24) we have 

^{z) = {z\^) = ^{z\j,m){j,m\^) 



= E 

m=-j 

= P^{z). 

The Majorana polynomial is therefore precisely the coherent state wavefunction. For each 
zero Zi of the coherent state wavefunction, its antipode — l/zj is the stereographic coordinate 



FIG. 1: A contour plot of y^\'>p{z, z)\ on the two-sphere for a typical spin state \ip). The location 
of its zeros is highlighted in red. The Majorana points of {ip) are the antipodes of these zeros. 

of the direction of one of the component spin-l/2's. These antipodes are known as Majorana 
points [16, 17]. 

In the following we sometimes find it convenient to use normalized coherent states 

and the corresponding non-holomorphic coherent-state wavefunctions 

iIj{z,z) =n{z\iP). 

In terms of these non-holomorphic wavefunctions the inner product becomes 

,,, , 2? + 1 /" dzAdz — , , . ^. 

The integral now involves only the usual area form on the two-sphere expressed in stereo- 
graphic coordinates. The normalized wavefunction iIj{z,z) is therefore the probability am- 
plitude for finding the spin-j pointing in the direction z. 

We can describe a spin-j state up to an overall phase by either specifying the zeros Zi of its 
spin coherent state wavefunction, or by specifying the Majorana points, i..e. the directions 
on the Riemann sphere of its component spin-1/2 factors. Because the spin directions are 
the antipodes Wi = —1/zi of the zeros, the normalized spin-1/2 state that points in the 
direction w has wavefunction 



This is just the wavefunction for the state lu?)^. Similarly we define 



\{Wi})j^ = \wi,W2,...,W2j) 



N 



= M ■ ■ ■ \w2j) (30) 

to be a normalized spin-j state with its 2j spin- 1/2 components pointing in the directions 
Wi. Here, Af is a Wi dependent normalization factor given by 

= = 7^ E + ^^^'^w)- (31) 

^ Tree i=i 

The corresponding wavefunction is 

2i 



AT 



{z\wi, W2,..., W2j)^ = (1 + l-^n-W"' JJ(1 + ZlDi). (32) 



1=1 



III. BERRY TRANSPORT IN THE MAJORANA REPRESENTATION 

Non-abelian Berry transport is a generalization of the usual Berry adiabatic transport to 
the case of degenerate energy eigenvalues [3]. Consider the adiabatic evolution of a state 

\^{t)) = an{t)\n,t) (33) 

that lies in a subspace V{t) spanned by a orthonormal frame \n,t),n — 1,...,N, of de- 
generate eigenstates of a Hamiltonian H(t). If H{t) varies slowly, the resulting adaibatic 
evolution leads to the coefficients being given by 

an{t) = Unm{t)am{0) cxp | -2 ^ E{t)dt^ , (34) 

where E{t) is the common eigenvalue, and 



Vexp< - / A 







(35) 



is a unitary matrix that generalizes the usual Berry phase. The symbol V indicates a path- 
ordered integral and the non-abelian Berry connection 

Anm = {n,t\d\m,t) = -Anm (36) 

is a skew-hermitian matrix- valued one-form. In the special case that degenerate subspace is 
one-dimensional it reduces to the abelian Berry connection. 
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We wish to find expressions for Anm in the Major ana parametrisation. The one dimen- 
sional (Abehan) case was obtained by Hannay [12]. He starts from 

= I^'^W) (^I^MO)) • • • ® k^(2j)) (37) 

where, from \wi) — Wi\'\) -\- \ \) we have d|wi) — diDilt), and 

^.n,n("f). (38) 



He obtains 



^ ^ ^ j^ E^ [nfc(l + WfcW^(fc)) Efc(wfcdw^(fc))/ (1 + WfcW^(fc)))] ^gg^ 

For J = 1/2, for example, we have 

f {w\d\w)\ l ivdw — wdw , , dw A dw , 

V {w\w) J 2 l + ww + ^ ' 

The expression for d^ is i times one-half the area form on the two-sphere, and so we find 

that familiar result that in a cyclic evolution the accumulated phase is one-half of the area 

swept out by the spin. 

For a general value of j consider the evolution of the E = m eigenstate of H = J ■ h 

as the unit vector n varies. In the Majorana parametrisation, j -|- m of the Wj's lie at the 

point w corresponding to the unit vector n and j — m of the Wi lie at the antipodal point 

w' — —1/w. Thus equation (39) becomes, 

1 , , wdw — wdw 1 , , w'dw' — w'dw , , , 

2 1 -I- WW 2 1 -I- w'w' 

For a cyclic evolution the accumulated phase is [(j + m) — (j — m)]/2 = m times the area 

swept out by the point w — again a familiar result [13]. 

Now consider the non-abelian extension of Hanay's formula in which Berry connection 

becomes a one- form A that takes its values in the Lie algebra u(n). We can evaluate its 

skew-hermitian matrix elements in the Majorana parametrisation by extending (39). After 

some work we find that the off diagonal elements of the connection are 

A (K}|d|K}) 

V({<}l{<})v/(K}IK}) 



^E. [n.(i + ^E. [n.(i + 
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and the diagonal elements are 

_ Im(K}|d|K}) 



({<}!{<}) 

(43) 





n.(i + «(.)) E.Kd<(,))/(i + «(,))^ 




n.(i+«(,))^ 





The expressions in equations (42) and (43) are, in general, rather complicated. Our interest, 
however, is in the evolution of Kramers-degenerate pairs. We therefore restrict ourselves to 
the time-reversal invariant hamiltonians, and to half-integer j for which the Kramers pair 
of states 1^) and T|^) are degenerate and mutually orthogonal. 
The time reversal operator on a general spin state is 

r = exp{-i7r J2}/C (44) 

where /C denotes complex conjugation. On a spin- 1/2 coherent state it acts as 

r\w) =w\- l/w) (45) 

and this action extends to the general Majorana-parametrisation tensor product state. 

Using the states l^'), Tj^') as basis, the Berry connection becomes a two- by-two traceless 
skew-hermitian matrix. 

The tracelessness of ^ is a consequence of time-reversal; it implies that the diagonal Abelian 
phases for a Kramers pair are complex-conjugates of each other. 
In the Majorana parametrisation 

^ ^ {{wi}\d\r{wi}) ^ Ett [nfc(w^(fc) - Wk) Efc(dwfc)/(w^(fc) - Wk)] 

' V{{M\{^i})V{r{wi}\T{w,}) E. [n.(i + wkw^ik))] 

(47) 

For eigenstates of H — J ■ h (47) simplifies considerably. Again, there are j + m Ma- 
jorana points w and j — m points at w' — —l/w. The numerator contains the factor 
rifc fc^j('"^7r(fc) ~ '^k) , which can only be non-zero for permutations that swap the w and w' 
except at most a single point that may map to itself (recall that 2j is odd). As a conse- 
quence the off-diagonal term is non-zero only for m = 1/2 states. When m ^ 1/2 the Berry 
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transport is purely Abelian, in the sense that there exists a fixed basis whose states are only 
multiplied by a phase for any loop in Hamiltonian space. 

Consider, for example rotations of n about the z axis. We have Wi — >■ Wiexpi(j) and 
dwi = iwid(j) and equation 47) becomes 

. Ep [Uki^^k) - Wk) Efc Wfc/(w^(fc) - Wk)] 

For a, m — 1/2 state with Majorana points at ty = tan(^/2) exp{i0} and w' — 
cot(^/2) exp{— i0} we get 

(n + l)!(n + l)! (1 + W)2"/^y'2n 
' ~ n!(n + l)! {1 + ww'f''+y{ww'y^'^^' 
1 

= (n+1)— -lyfid^) 

= exp [i(2n +1)0] ^^^^^ sin(^) (id0) (49) 

where 2n + 1 = is the total number of Majorana points. We may redefine the basis states 
so as to absorb the phase factor so that 

^1 = sin(e)(id0) = ^ (j + sin(e)(id0), (50) 

and this expression agrees with a result in [5] . 



IV. KATO'S EQUATION 

The Berry-connection one-form allows us to adiabtically evolve states, but, because it 
involves continually making a choice of basis, it is not always the most convenient method 
for numerical computation. An alternative tool is provided by the operator Pif) that projects 
us into our space V{t) of eigenstates of Hif) with a common eigenvalue E{t). This subspace 
evolves with the parameters in i?(i), and the Berry-transport condition [3] is 

n^)|w^)) = 0. (51) 

Intuitively, the Berry condition says that we keep projecting the state down into the evolving 
degenerate subspace. This process is completely equivalent to transport via Anm- 
Using the projection operators we can approximate the adibatic evolution as 

m)) = P{t)P{t - e)P{t - 26) • • • P(0)|V^(0)). (52) 
13 



For finite e, this sequence of successive projections is not norm-preserving, but in the hmit 
e — > it becomes so. 
We now observe that 

-^{m+e))-m))} = l{p{t+e)-i}p{t)p{t-e)---pmm) 

= ^{P{t + e)-P{t)}P{t)P{t-e)---P{Q)\m- (53) 

Here we have used the projection-operator property P{t)'^ = P{t)- On taking the limit 
e ^ 0, we are lead to evolution equation 

^^\m) = pmit))- (54) 

Let us confirm that this evolution satisfies equation (51): from P'^ — P we have PP + PP — 

P, and so 



d_ 

= {pp-Pm 



P-M)t = PPit)\ym 



= {P-PM 

= 0. (55) 

In passing from the second line to the third we have used PiV') = IV') for IV') iii the degenerate 
subspace V{t). 

Although equation (54) is very simple, it has a formal drawback in that it does not provide 
a manifestly unitary state evolution. The operator P{t) is hermitian, and so therefore is 
P{t). For unitary evolution we would expect the operator on the right hand side of (54) to 
be skew hermitian. We can, however, equivalently write (54) as 

^^m,^[p{t),p{t)]m))- (56) 

We have already seen that the the added term PP vanishes when acting on states in the 
degenerate eigenspace. Now the commutator of two hermitian operators on the right-hand 
side is skew hermitian, so this version of the evolution equation does provide a manifestly 
unitary evolution. Equation (56) governing the non-abelian Berry transport is known as 
Kato's equation [15] after Tosio Kato who obtained it in 1950. 
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Now our interest is in the evolution of the spin in a molecule that is rotating in space. 
For a fixed orientation of the molecule, the spin Hamiltonian is a polynomial 

2j 3 

H(^)^Yl Yl (^hi2-imJilJi2---Jim (57) 

m=0 ii,...im=l 

in the space-fixed components Jj of the spin angular momentum operator. The effect of a 
real space rotation R on the molecule is to is to change Jj — > R^j^Jj in this polynomial while 
keeping the coefficients ajjj2...j^ fixed. Now for each it! there is a unitary operator U[R\ that 
implements the rotation in the 2j + 1 dimensional Hilbert space: 

R-f'jj^U[R]JiU-'[R], (58) 

so as the molecule rotates we can write 

H(t) = U[R{t)]H{0)U[R{t)]-\ (59) 

where U (t) — U[R{t)]. (This U (t) is not to be confused with the adiabatic evolution operator 
Unm{t) in equation (35).) The same operator evolves the projection 

P{t) = U{t)P{Q)U{t)-\ (60) 

and we can work in a frame rotating with the molecule by setting \4'{t)) = U{t)\(f){t)) . In 
this frame Kato's equation becomes 

The evolution given by (61) now all takes place within V{0) and = G{t)\(j)o), where 
the propagator G{t) is be the time-ordered exponential 

G{t) = Vexp(^- P{0)K{t)P{Q)ds 

= lim {P(0)(1 - eK(t))P{0)(l - eK(t - e))P(0) • • • P(0)(1 - eX(0))P(0)} , (62) 

where K = U'^U . From (62) we can see that G{t)P{Q)G{t)-^ = P(0). 



V. ACTION PRINCIPLE AND EQUATIONS OF MOTION 

The projection P(0) in (62) acts as the identity operator on the subspace V{Qi), and is 
effectively the quantum evolution of a state restricted to V{Qi). We can therefore represent 
(62) by means of a Feynman path integral over trajectories in some suitable classical system. 
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As in the previous section we assume time-reversal invariance and half-integer spin, so 
that by Kramers' theorem each energy level is doubly-degenerate and spanned by |0(O)) and 
T|0(O)) (which we denote here as |T0(O))). In such a space the initial projection operator 
P(0) is 

p{o) = \m){m\ + \rm){rmi m 

Since V{0) is two-dimensional, the natural path integral is the spin-coherent-state path 
integral for SU(2). 

To obtain the path integral we let A G CP^ be the stereographic complex coordinate of 
a point on the Riemann sphere and set 

|A)^ = (l + AA)-i/2(|r0(O)) + A|0(O))) (64) 
(A|^ = (l + AA)-V2((r0(O)| + A(0(O)|). (65) 

These are analogous to the spin-1/2 coherent states in which |T0(O)) and \4>{0)) are playing 
the roles of the spin down and spin up states respectively. 

We now write the projection operator P(0) as the (overcomplete) resolution 

^/^\ 2 f 2dAAdA ,,,,,, 

^(°' = sy„.i(TTW"><"^' '''' 

and with it construct, by time-slicing with P(0), what is essentially a spin-1/2 "coherent" 
state path- integral [11] containing the action 

S[X, A]] = £ dt (^l^f^ - A(A, A) j , (67) 

with Ao = {X\K\X). 

We might attempt to evaluate this path integral by some technique. However, when the 
Hamiltonian term is linear in the Lie algebra generators, as it is in this case, then we know 
that the classical equations of motion give the exact quantum evolution (See for example 
[18, 19]). 

These equations are 
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Now, from the construction of |A), we can write 



ii + xx)t{x\x\x)^-{x\xr\x) (70) 

for any operator X. In particular we can rewrite (69) as 

^ = -(1 + AA)A (71) 
^ = -(1 + XX)A^ (72) 

where ^i(A, A) = {X\KT\X). We now know the time evolution of X{t), and hence the time 
evolution of |0(i)). 

To see how changes in A and A translate into changes in the location of the zeros Zj, or 
equivalently the Majorana points Wi = —l/zi we observe that 

|A)^ = e'1W})^ (73) 

for some phase 7, and so 

ilj{z,x) = e'^^^-\l + \z\Y'Yli'i- + zwi) 

i=l 

^Ql[{l + zw,), (74) 

i=l 

where Q = e^'^N'~^{l + |zp)~-'. The last form reminds us that the roots of the wavefunction 
i(j{z, A) are at Zi = — l/wj. Assume for now that the zeros are all non-degenerate — that is, 
Zi 7^ Zj if i ^ j. Consider a root Zk{X) of {z\X), and recall that the zeros of {z\X) depend 
only on A while those of {X\z) depend only on A. Let A — > A + 5A and expand ^^(A + 5X) to 
first order as 

z,CX + SX)^z,{X) + ^5X, (75) 

dX 

For Zk{X + SX) to remain a root, we must have 

ip{zk + -^5X,X + 5X) = Q. (76) 
dX 

Expanding (76) to first order and using ip^Zk-, A) = 0, we find [6, 20] 

(77) 



dzk dif) f dip 
'dX ^ ~dX [~d^ 
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The denominator of (77) is obtained by straightforward differentiation of (74) to be 



dz 



= -Qwl-''ll{wk-Wi), (78) 



and the numerator is 

= -(l + lXiye-^'-'Qwf^Ylil + WiWk). (79) 



dip 
dX 



Combining equations (76), (78), (79), and using that 

dzk d f 1 \ 1 dwk 



dX dX \ Wk J wl dX 
we arrive at differential equations for the Majorana points Wk in terms of A, 



(80) 



(l + AA)—— = e + , (81) 

dX Wk - Wi 

(1 + AA)— ^ = e ^{l + WkWk)\\— —. (82) 

dX Wk - Wi 

Now 

dw dw dX dw , ,,,9s . /„„s 

where we have used (72), and 

A = n{X\U-^UT\X)^ = e-^''^^{{w}\U-^Ur\{w})^. (84) 
The phase 7 cancels leaving us with 

dWk .r 1- , _ -p-r 1 + WiWjfc 
-— = JV({W}|[/ r|{w})^(l + WfcWfe) I I , 

at " Wk — Wi 

ij^k 

^ = ^r(W|C/-^?7r|W)^(l + «) n (85) 
at f Wk — Wi 

The equations (85) are quite intricate. To gain some insight, consider some of their 
properties. We begin by locating possible fixed points of (85). If Wk is to be zero there 
are two possibilities: either the the product rii^ifc(l + 'WiWk){wk — Wi)~^ must vanish, or 
the factor n{{w}\U~^UT\{w})j^ must vanish. We focus on the first possibility. For this 
to occur, one of the Majorana points Wi must become antipodal to Wk — that is there is a 
Wk'it) such that 

Wk'{t) = (86) 
18 



At that instant, Wk — 0. But observe that this also hold true under k — )■ k', and Wk — as 
well. Consequently, both Wk and Wk' remain static for all t. As such, antipodal Majorana 
points are always fixed. 

We can understand this by noting that antipodal Majorana points are actually Majorana 
points that are common to both and \T4>t)- Suppose the Majorana points of \4>t) are 
given by 

ai,-l/ai,...,a^,-l/a,.,^i,...,^p (2r+p = 2j) (87) 

where the Oj's come in antipodal pairs. Since T takes the Majorana points to their antipodes, 
the Oj's are preserved; in other words, the Majorana points of \T4>) are 

ai, -1/ai, . . . , a^, -1/a^, -1/6, • • • , -l/Cp- (§8) 

Note that 10^) is a linear combination of |0o) and \T(/)o), when 10^) and |T0o) have any 
Majorana points in common, then will also share these Majorana points. 

To summarize: if a pair of Majorana points start out being antipodal, then they remain 
fixed at all times. In future we will factor out these fixed points, and refer to the non-fixed 
Majorana points as "free" or "dynamical" . 

Since the fixed Majorana points do not participate in the dynamics, the requirement in 
equation (85) of non-degenerate Majorana points can be relaxed to only requiring that none 
of the free Majorana points ever coincide. This includes states whose fixed antipodal pairs 
occure more than once (for example, the \j, 1/2), \k, —1/2) states). 

Writing ^ for the free Majorana points and Oj, — 1/aj for the fixed antipodal pairs as in 
(87), we find that 

^ = N{{a},{m-'Ur\{a},{0)^{l + ^k^k)v n 

^ = N{{a},{0\U-'UT\{a},{0)^{l + ^k^k)v JJ (89) 

Here rj — YVi=ii~^i/'^i) a-dependent overall phase, and denotes a state in 

the Majorana representation with fixed Majorana points a, — 1/a, and free Majorana points 
6 Observe that (89) is well behaved even if coincides with one of the a's, but becomes 
singular if any ^ approaches another free ^k' ■ 
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An especially interesting case occurs when only one Majorana point is free and the rest 
compose antipodal pairs. In this case the equations of motion reduce to 

^ = N{{a},c\u-'ifr\{a},c)Ai + COv (90) 



§ = N{{a},m-'Ur\{a},ONi^ + ik^k)v (91) 
What is interesting here is that (91) can be derived from an action 

S[C,C] = j\s (^l^i^-N{{a},m-'U\{a},ON^ ■ (92) 
The semiclassical equations of motion obtained from this action are 

§ = -(l + ay^iv({a},^|^/-^^7|{a},0^. 

{i + ^0'§iN{{a},m-'U\{a},ON, (93) 

At first glance equations (91) and (93) look different. That they are, equivalent however, 
can be verified by using equation (70) 

VI. ILLUSTRATION 

As an application of our formalism, consider the family of molecular magnets Mn4. The 
effective Hamiltonian for the j — 9/2 molecular spin in the absence of an external field can 
be written as 

n = + kyJy + kzJz, (94) 

together with small quartic terms that we will ignore. (The small field-coupling between 
the spins on different molecules will also produce a term linear in J — and therefore time- 
reversal symmetry breaking — which may be significant. We are, however, also going to 
ignore this for the moment.) Following [22], we adopt the values — —ky — 0.0243K, 
K = -0.59K. 

First, we look at the configurations of the zeros for the eigenstate wavefunctions of %, 
for which there are a total of 5 doubly-degenerate levels when j = 9/2. Due to the relatively 
large quadrupole moment kz of the z axis, the Hamiltonian is close to that of J|, and as a 
result the zeros of the wavefunctions will want to cluster about the poles (see figure 2). But 
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FIG. 2: The zeros of the eigenstate wavefunctions for 7i in absence of an apphed field, for which 
the expectation (J) is largest about the z-axis. Blue and red denote the zeros of the eigenstate and 
its Kramers pair, repsectively. In order from left to right and top to bottom, we have (1) ground 
state, (2) first excited state, (3) second excited state, (4) third excited state, (5) highest state. In 
the background we have plotted the energy contours of corresponding semiclassical Hamiltonian. 
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unlike zeros of {z\j, m) which simply condense at the poles, the small and Jy terms will 
cause some repulsion between them. In the background we also display the contours of the 
"classical" potential landscape defined by the expectation h{z, z) =^ m{AH\z) j^. 

Due to the double degeneracy however, there is a freedom in the choice of the basis of the 
eigenspace. One particularly illuminating selection is made by putting in a small transverse 
magnetic field in the x direction and relaxing the field strength to zero. In figure 3 we see 
that for the ground states the zeros tend to distribute themselves along flow hues Vh{z, z) 
in a manner such that the two ground states have all their zeros distributed about regions 
of maximum energy, while for the first excited states a pair of zeros for each state have 
relocated themselves close to the minima, for the second excited state an additional pair, 
etc. In addition, for the first, second, and third excited states we find that the zeros about 
the poles of one state very closely overlap with that of its Kramers pair, and this can have an 
effect when considering the dynamics their dynamics under adiabatic evolution. We remark 
however that for this choice of states the magnetization (J) becomes significantly reduced, 
and may consequently be very difficult to measure. 

Now we consider adiabatic holonomy arising from rotations about a fixed axis fi = 
sin?7cos7x + sin 77 sin 7 y + cost] 2 at an angular velocity of l-n jT . For such rotations we 
are guaranteed to have the state return to itself (in the corotating frame) after some time 
T, which need not coincide in general with the period of rotation T; the former is gotten by 
decomposing PU~VP into sigma matrices Xiai in the eigenstate basis, for which the period 
is then found as 27r/\x\. Fixed-axis rotations have been considered before by [5], but their 
attention is restricted to eigenstates of isotropic Hamiltonians of the form n • J. 

The motion of the zeros can be quite complicated. The greatest intricacy occurs for higher 
levels where the Kramers pairs communicate more strongly. As an illustrative example we 
plot in figure 4 for one full cycle r, and for various axis orientations, the trajectories in 
the corotating frame of the highest- level wavefunction zeros (polarized initially along x). 
Observe that trajectories of nearby axes can exhibit qualitatively different behavior — in 
particular closed cycles of individual zeros can combine into larger cycles in which the zeros 
permute locations. These trajectories have been calculated numerically by both using the 
projection method in Kato's equation, as well as by solving the equations of motion for the 
Majorana points, and are found to agree well. (The intricate motion is best appreciated by 
viewing animations. We have made these available at [21]) 
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FIG. 3: The zeros of the eigenstate wavefunctions for 7i in absence of an apphed field, for which 
the initial states are prepared by polarizing along the x. Note that in the first, second, and third 
excited states the zeros about the poles actually very closely overlap — for example the blue dot 
in the north pole for the first excited state is actually a pair of blue and red dots. 



23 



FIG. 4: Trajectories of the zeros, in the corotating frame, for the highest energy level. Lighter 
blue corresponds to earlier in time, and darker to later. In the left figure the axis was chosen along 
rj = 0.0345-7r, 7 = 0; and in the right figure the axis was along rj = 0.03557r, 7 = 0. A small change 
in zenith angle has completely changed the topology of the trajectories. 

Now consider the time- dependence of the magnetization (J). By spinning the molecule at 
some frequency 1/T about the chosen axis, the magnetic moment may oscillate at a different 
frequency 1/r, which can be measured by placing pickup coils near the sample. The exact 
form of the moment oscillation depends on the energy level as well as the choice of initial 
state within the level. 

In figure 5 we look at a sample initially in the ground state, polarized along the z direction. 
Such a state has a large initial magnetization along z and zero magnetization along the 
transverse directions. We then rotate the sample about x at a frequency of 2tt/T, where T 
is large compared the timescale of the gap — in this case «i 8 K, or ^ 6 x 10~^^ seconds. 
We then calculate (J) in the body frame of the molecule, and we find that while the (J^) 
and (Jy) components do not vary much (being virtually zero), the (J^) component exhibits 
appreciable oscillations on the timescale of r 165000T. What one should see in the lab 
frame then are rapid oscillations of (J^) on the order of T, with an amplitude modulation 
with period r. 

Finally, in figure 6 we look at a state initially polarized along the x direction and rotated 
about z. The motivation for this is to look at states whose magnetization along the axis 
of rotation is suppressed. We find that as the energy level becomes higher, the component 
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Expectation of Angular Momentum in Body Frame 
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FIG. 5: Plot of (J(t)) for the ground state polarized along z, with {J^) in red (not visible), {Jy) in 
green, and (J^) in blue. 

{Jz) of the initial state becomes rapidly smaller , while (J^) grows. We expect that at 
higher levels non-Abelian effects become stronger. Indeed we find that in the body frame, 
oscillations of the ( J^,) and {Jy) components become more pronounced, which may be easier 
to detect. It would therefore be interesting if one can select an initial energy level, perhaps 
by applying a suitable a/c pulse to excite the molecule and then selecting the starting state 
within that level by applying a small d/c field which is gradually reduced to zero. 

VII. DISCUSSION 

We have derived the equations of motions of the Majorana points for states undergoing 
non-Abelian adiabatic evolution. We focussed primarily on the case when the degenerate 
subspace consists of a Karmers pair of eigenstates for a time-reversal invariant, half- integer 
spin system, and on time-dependence that arises from simply rotating the host molecule 
about a fixed axis. We have seen that even when the effective Hamiltonian is quite simple 
(quadratic anisotropy only) the resulting dance of the Majorana points can be quite intri- 
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FIG. 6: Plot of (J(t)) for the third (top) and fourth (bottom) excited states polarized along x, 
with {Jx) in red, {Jy) in green, and (J^) in blue. 
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cate. When only one Majorana point is able to move its equation of motion is manifestly 
Hamiltonian. When many Majorana points are in motion the system is still Hamiltonian, 
but the associated symplectic form will be complicated. 

Non-Abelian evolution means that the motion of the period of the Hamiltonian (the 
rotation period) and the period of the spin direction will in general be incommensurate. 
The resulting difference between the input rotation frequency and the output magnetic 
moment frequency can in principle be measured. A confounding effect, however, is provided 
by the very moment that we are trying to measure. If it is too large, the resulting magnetic 
field will destroy the degeneracy of the Kramers pair and will wipe out the non-Abelian 
dynamics. If it is too small there will be no signal indicating the spin direction. The 
splitting of the degeneracy means that the two not-quite degenerate states will acquire a 
steadily accumulating relative phase that will ultimately overwhelm the non-Abelian affect. 
To avoid this problem we need to have a system in which there is a substantial difference in 
scale between the small splitting of the near-degenerate states and the gap to the first excited 
state of the time-reversal invariant system. We can then spin the system fast enough that 
the small splitting can be neglected, while the evolution within the Kramers pair eigenspace 
is still adiabatic. 

VIII. ACKNOWLEDGEMENTS 

The work of YL was supported by grant NSF PHY08-55323AR, and that of MS and AR 
by grant NSF DMR 09-03291. We thank Anupam Garg and FeiFei Li for useful discussions. 



[1] W. Wernsdorfer, L. Bogani, Nature Materials 7, 179 - 186 (2008). 

[2] W. Wernsdorfer and R. Sessoli, Science 284 133-135 (1999). 

[3] F. Wilczek, A. Zee, Phys. Rev. Lett. 52 2111-2114 (1984). 

[4] J. Moody, A. Shapcrc, F. Wilczek, Phys. Rev. Lett. 56 893-897 (1986). 

[5] C. A. Mead., Phys. Rev. Lett. 59 161-164 (1987). 

[6] P. Leboeuf, J. Phys. A: Math. Gen. 24 4575-4586 (1991). 

[7] E. Majorana, II Nuovo Cimento, 9 43-50 (1932). 

[8] J. Schwinger, Phys. Rev. 51 648-651 (1937). 

27 



[9] 1. Rabi, Phys. Rev. 51 652-654 (1937). 

[10] F. Bloch, I. Rabi, Rev. Mod. Phys. 17 237-244 (1945). 

[11] M. Stone, K-S. Park, A. Garg, J. Math. Phys. 41 8025-8049 (2000). 

[12] J. H. Hannay, J. Phys. A31 L53-59 (1998). 

[13] M. V. Berry Proc. Roy. Soc. Lond. A 392 45-57 (1984). 

[14] J. E. Avron, L. Sadun, J. Segert, B. Simon, Comm. Math. Phys. 124 595-627 (1989). 

[15] T . Kato. J. Phys. Soc. Japan, 5 435-439 (1950). 

[16] M. Aulbach, D. Markham, M. Murao, Proceedings of TQC 2010, Lecture Notes in Computer 

Science 6519 141-158 (2011). 

[17] J. Zimba, European J. Theor. Phys. 3 143-156 (2006). 

[18] G. Giavarini, E. Onofri J. Math. Phys. 30 659-663 (1988). 

[19] M. Stone, Int. J .Mod. Phys. B15 2107-2114 (2001). 

[20] D. Elhnas, V. Kovanis, Phys. Rev. A51, 4230-4239 (1995). 

[21] http: / / mstone.physics.ilhnois.edu/ 

[22] F. L. Mettes, G. Aromi, F. Luis, M. Evangehsti, G. Christou, D. Henderickson and L. J. de 
Jongh, Polyhedron, 2011, 20, 1459-1463. 



28 



